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This paper presents a theoretical treatment of the roundoff noise problem 
for the special case of cascade realizations of Finite Impulse Response 
(FIR) digital filters. 1 Explicit relations for evaluating roundoff noise with 
the usual assumption of uncor related samples are presented. Useful scaling 
methods are stated and classified as to conditions when these methods are 
optimum. Important differences between use of these scaling procedures for 
Infinite Impulse Response (IIR) filters and FIR filters are pointed out. 
Finally, useful properties of linear phase FIR filters are discussed. 

I. INTRODUCTION 

In recent years, many techniques have been developed for the 
design of Finite Impulse Response (FIR) digital filters. 2-8 It is now- 
possible to readily design filters with arbitrary frequency or time 
response characteristics using the windowing, 23 frequency sampling, 4 
or optimal design 5-8 methods. While both the windowing and fre- 
quency sampling techniques yield suboptimal filters, they are useful 
because of their simplicity and ease of design. The optimal design 
technique is of special importance because the filters it generates can 
be proved to be optimum in a certain sense, 7 and because efficient 
algorithms exist for its implementation. 7 ' 8 

As a result of these important developments, the FIR type of 
digital filter is becoming increasingly attractive as an alternative to the 
IIR (Infinite Impulse Response) type of filter for practical applications. 
A major advantage of FIR filters over IIR filters is that an FIR filter 



* This paper is based on a thesis' submitted in partial fulfillment of the require- 
ments for the degrees of Bachelor of Science and Master of Science in the Department 
of Electrical Engineering at the Massachusetts Institute of Technology in September 
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can have an exactly linear phase response while approximating an 
arbitrary magnitude frequency response. But even without considering 
this important advantage, current research 9 is revealing that in certain 
cases FIR niters are competitive with IIR filters in terms of speed and 
cost. Thus the implementation of FIR filters using finite-precision 
arithmetic is becoming an important area for research. 

Up to the present, little is known as to how different types of FIR 
filter realizations behave with respect to quantization effects. Since 
hardware, specifically for the purpose of realizing FIR filters, has 
already been built for experimentation by various research groups, 10 - 11 
it is important to obtain more knowledge to guide the implementation 
phase of FIR digital filter design. The purpose of this paper is to 
present a theoretical treatment of several problems associated with 
implementations of these filters. 

II. PRELIMINARY REMARKS 

The effects that quantization has on an IIR filter can be classified 
into three basic categories : 

(i) Quantization of the values of samples derived from a con- 
tinuous input waveform causes inaccuracies in the representa- 
tion of the waveform (A-D noise) . 
(it) Finite-precision representation of the infinite-precision filter 
coefficients alters the frequency response characteristics of the 
filter (coefficient accuracy problem). 
(Hi) Finite-precision arithmetic causes inaccuracies in the filter 
output (roundoff noise) which, together with the finite dynamic 
range of the filter, limit the signal-to-noise ratio attainable. 
Also, finite-precision arithmetic can lead to limit cycles where 
the output samples are generally highly correlated. 

These same quantization effects also occur in finite wordlength FIR 
filters with the important exception that limit cycles cannot occur in 
nonrecursive realizations of FIR filters. In this paper only the third 
type of quantization effect, viz., roundoff noise, will be discussed. 
Furthermore, fixed-point arithmetic with rounding will be assumed. 

Except for the first category above (A-D noise), all quantization 
effects depend in degree and character on the type of structure used to 
implement a filter. There are three well-known structures in which an 
FIR transfer function can be realized. They are the direct form, the 
cascade form, and the frequency-sampling structure. 12 Other less 
well-known structures based on polynomial interpolation formulas 
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Fig. 1 — Cascade-form filter section. 



are also possible; these include the Lagrange, Newton, Hermite, and 
Taylor structures. 13 However, it is as yet unclear under what circum- 
stances, if any, these other structures may be more advantageous than 
the well-known structures. 

Only the cascade structure will be discussed in this paper. A second- 
order filter section, as shown in Fig. 1, will be used as the basic building 
block for the cascade structure. Although several minor variations to 
this configuration for the filter sections are possible, 1 the results pre- 
sented here are sufficiently general so that they can be readily applied 
to other configurations as well. 

Aside from section configuration, the prime issues that must be con- 
fronted in the realization of filters in cascade form are scaling and 
section ordering. Proper scaling must be performed on a filter in 
cascade form in order that full use of the dynamic range of each section 
can be made while avoiding error-producing overflows. By proper 
scaling, the signal-to-noise ratio of a filter can be maximized for a 
given quantization step size and section ordering. 

Proper ordering must also be determined for a filter in cascade form 
if the filter is to be useful at all, since the noise output of a cascade 
filter can depend dramatically on the way it is ordered. For example, 
Schussler 13 showed a 32nd-order FIR filter which, ordered two different 
ways, yields output noise variances that differ by a ratio on the order 
of 10 8 . The problem of section ordering for cascade FIR filters has been 
investigated in depth, 114 and the results show that for higher-order 
filters, the variation of output noise variance across all orderings is 
much greater than 10 8 . 

Jackson 15 has formulated the roundoff noise problem for a general 
digital filter and has proposed an approach to the scaling of filters to 
satisfy dynamic range constraints. Most of his results can be specialized 
to the case of FIR filters by assuming a constant polynomial in the 
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denominator of the transfer function. However, the different perspec- 
tive obtained by studying the FIR case separately affords many 
additional insights. 

In this paper, formulas for the evaluation of roundoff noise variances 
in the FIR cascade structure are presented. Also, specific scaling 
methods for FIR filters are denned, two of which can be proved to be 
optimum for two classes of input signals. Finally, certain properties 
of the linear phase cascade structure useful in the study of section 
ordering are stated and can be proved. However, for reasons of space, 
no proofs are included in this paper. They can be found in Ref. 1. 

III. DEFINITIONS 

The general transfer function for an JV-point FIR filter can be 
written in the form 

H(z) = £ Hk)z~ k (1) 



k=0 



where the real-valued sequence [h(k), k = 0, • • •, N - 1} is the im- 
pulse response of the filter. Alternatively, H(z) can be expressed in the 
factored form 



-V.v 



H(z) = II (&o< + buz' 1 + & 2 rfT 8 ) (2) 

where bji, j = 0, 1, 2, i = 1, •• ■, N a are real numbers and N„ the 
number of factors, is defined as 

cN - 1 
N, = • 



2 


N odd 


N 
.2 


N ever 



and b 2N> = if iV is even. 

A linear phase filter is defined to be a filter, the transfer function H (z) 
of which is expressible in the form 

ff(«)|,_.y« = H(e^) = ±\H(e*")\e-* a ° (3) 

where a is a real positive constant with the physical significance of 
delay in number of samples. The factor ± is necessary since H(e ju ) 
actually is of the form 
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where H*(e jw ) is a real function taking on both positive and negative 
values. It is useful to define a mirror-image polynomial (MIP) of 
degree N to be a polynomial of the form £jlo a k z k , the coefficients of 
which satisfy the relation 

a* = a,N-k ^ k ^ N. 

Necessary and sufficient conditions on H(z) such that a filter with 
transfer function H(z) has an exactly linear phase response can then 
be stated as follows : 

Theorem 1: H(z) can be expressed in the form (3) if and only if one of 
the following equivalent conditions hold: 

(i) h(k) = h(N - 1 - k), ^ k ^ N - 1. 
(it) If Zi is a zero of H(z), then z~ l is also a zero of H(z). Also, if 

2, = + 1 is a zero of H(z), then it occurs in even multiplicity. 
(Hi) Suppose Zi is a zero of the ith factor in (2). Let S = [i: Zi is real] 
and Q= [i:i$S\. Thenf(z) = U<es (&« + b u z~ l + b u z~ 2 ) is 
a mirror-image polynomial in z _1 , and for all i £ Q, either 
boi = &2. or there exists j ^ i, j £ Q, such that 



bo, _ bu b 2 i 
&2> bu b j 



(4) 



Furthermore, the following is a sufficient condition for H(z) to be expres- 
sible in the form (3) : 

(iv) In (2), for 1 ^ i ^ N„ either b u = and b 0i = b u , or b 0i = b 2i , 
or there exists j j* i, 1 ^ j : ^ N„ such that 

boi b u b 2 i 

b 2} bu b j 

In all cases the value of a is a = (N — l)/2. 

It should be pointed out that a section with 6 , = 621 is necessarily 
one which synthesizes either two complex conjugate zeros on the 
unit circle, or two reciprocal zeros on the real axis, or two identical 
zeros at +1 or — 1. Furthermore, two sections which satisfy (4) are 
precisely those sections which synthesize reciprocal zeros (i.e., if 2, 
is a zero of one section, then zf 1 is a zero of the other section). Thus, 
taking (2) as the basis for the FIR cascade form, condition (iv) of 
Theorem 1 provides a way to assign zeros to individual sections of 
the cascade structure so that linear phase is guaranteed independent 
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of scaling or ordering. Hence, in this paper the following convention 
of zeros assignment for linear phase filters will be adopted : complex 
zeros are grouped by conjugate pairs, real zeros that are reciprocals of 
each other are paired together, while doubled or higher multiplicity 
zeros are grouped by pairs of the same kind. In this way the only zero 
that can occur by itself in a section is z = - 1 (since by Theorem 1, 
2 = + 1 is not allowed as a zero of odd multiplicity). 

The definition (3) of a linear phase filter requires the filter to have 
both constant group delay and constant phase delay. However, if 
only constant group delay is desired, a second type of "linear phase" 
filter can be defined in which the phase of H(e s ") is a piecewise linear 
function of w, i.e., 

H(e*") = ±|#(e''»)|e' , «»-«» ) . (5) 

It can be shown 1 that with the constraint (1) on the form of H(z), the 
only possible solutions for e [-*, x] is = ± (fcir/2), k = 0, 1, 2. 
If = 0, ±tt, (5) reduces to (3). Thus the only new cases added are 
when (8 = db tr/2. These cases arise exactly when Zi = + 1 occurs as 
a zero of H(z) in odd multiplicity, or, equivalently, when \h(k)\ 
satisfies 

h(k) = - h(N - 1 - k) £ k g N - 1. 

Filters of this special type are useful in the design of wideband differ- 
entiators. 16 However, this type of filter will not be considered in this 
paper and the term "linear phase filter" will be restricted to refer to 
those filters satisfying (3). 

IV. THEORY OF FIR CASCADE STRUCTURES 

4.1 Roundoff Noise in the Cascade Structure 

The analysis of roundoff noise in this paper is based on the usual 
model used for such analyses in digital filters. 15 ' 17 ' 18 In particular, each 
multiplier in a filter is modeled as an infinite-precision multiplier 
followed by a summation node where roundoff noise is added to the 
product so that the overall result equals some quantized level. Each 
noise sample is modeled as a random variable with uniform prob- 
ability density on the interval (-Q/2, Q/2) and zero density elsewhere, 
where Q is the quantization step size. Thus each sample is a zero-mean 
random variable with a variance of Q 2 /12. 

Furthermore, the following assumptions are made : 

(t) Any two different samples from the same noise source are 
uncorrelated. 
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(ii) Any two different noise sources (i.e., associated with different 

multipliers), regarded as random processes, are uncorrelated. 
(Hi) Each noise source is uncorrelated with the input signal. 

Thus each noise source is modeled as a discrete stationary white 
random process with a uniform power density spectrum of magnitude 

QV12. 

Applying this model to the filter section shown in Fig. 1, the addition 
of a noise source to the output of any multiplier is seen to be equiv- 
alent to adding a noise source to the output of the section. Therefore, 
to model a section of a cascade filter, fc, noise sources are added to the 
output of the section, where k t is the number of multipliers with non- 
integer coefficients in the section. Or, equivalently, by assumption (ii) 
above, one noise source of variance A;,-(Q 2 /12) can be added instead. 

For the configuration shown in Fig. 1, fc, is in general equal to 3. 
However, when b > = b 2i , the signals of the two branches feeding the 
multipliers with coefficients b 0i and b 2 . can first be summed before 
being multiplied by the common coefficient, thus reducing A;,- to 2. 
Furthermore, by a sacrifice in speed (assuming serial arithmetic), 
it is possible, as demonstrated by practical hardware, 10 to reduce 
A;, to 1 for all i by summing all products in each section before per- 
forming rounding. It is of interest to point out that the same can 
be done in the direct form, resulting in effectively only one noise 
source of variance Q 2 /12 feeding into the output of the filter. 

Before proceeding further, some notations need to be developed. 
Let Hi(z) denote the transfer function of the fth section of a filter 
H(z), i.e., 

H(z) = if H&) (6) 

where 

Hi(z) = b 0l + b u z-' + b u z- 2 . (7) 

As a convention, filter sections will be numbered in increasing numbers 
according to increasing distance from the filter input (i.e., the section 
at the input is called the 1st section). 

Furthermore, define 

II H,(z) < i < N, - 1 
G,(*)= ,-T+i " " (8) 

I 1 i = N, 

and let \gi(k)\ be the impulse response of Gi(z), i.e., 

Gi(z) = £ 9i(k)z-». (9) 

k 
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Fig. 2 — Equivalent models for a filter in cascade form. 

Then a cascade filter can be modeled as in Fig. 2a or equivalently 
as in Fig. 2b, where {e,-(n)| is the noise source for the ith section. 
Letting [Ei(n)\ denote the noise sequence at the filter output due to 
the zth noise source alone gives 

EM = ZgtikMn- fe). (10) 

k 

By the stationarity of {e,(n)}, the variance of Et(n) is independent 
of n ; hence denoting this variance by a? , assumption (i) above leads 
to the relation 



<n = L gi{k)ei{n - k) 2 



12 k 
Now the total noise output is given by 

E(n) = E Ei(n) = E E g<(fc)e,-(n - fc). 



(11) 



(12) 
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Therefore by assumptions (i) and (it) 



N. 



ff 2 _ E 2 {n) = £ „ im (13) 



i = l 



4.2 Methods of Scaling to Meet Dynamic Range Constraints 

A practical digital filter, necessarily implemented as a physical 
device, must have a finite dynamic range. Especially when fixed-point 
arithmetic is employed, this dynamic range sets a practical limit to 
the maximum range of signal levels representable in a filter and acts 
to constrain the signal-to-noise ratio attainable. 

In some filter structures such as the direct form, given the filter 
transfer function, the designer has no control over the relative signal 
levels at points within the filter. Only the gain of the overall filter 
can be varied. However, in a cascade realization with N s sections 
there are N, — 1 degrees of freedom available in addition to the 
overall filter gain and the ordering of sections. 

To see this, a factorization for H(z) is defined which is unique up 
to ordering of factors, in the form 

h(z) = n aw 

Hi{z) = a 0i + ausr 1 + a zi zr* (14) 

where {a,->| satisfies 

2 

oo,- e=0, £ M =1 i = 1, ••• ,N.. (15) 

;=0 

Then the transfer function for the ith section in a cascade realization 
can be written as 

Hi(z) = Sifiiiz) (16) 

where Si is an arbitrary constant, subject only to the constraint that 

II Si = 0. (17) 

Thus given /3, N„ — 1 of the Si 6 can be chosen at will. 

Any rule for assigning values to [St] will be referred to as a scaling 
method. Obviously, some scaling method must be employed in the 
design of a cascade filter whether or not one is concerned with dynamic 
range constraints, since numerical values must be assigned to the Si's. 
When dynamic range is an issue, the constraints it imposes can be 
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met in some best manner by choosing the proper scaling method. 
In this paper, filters, designed so that no arithmetic overflow in them 
can cause distortion in the filter output, will be studied. Therefore, the 
investigation of scaling methods will be restricted to those methods 
which guarantee that for a given class of input signals no distortion- 
causing overflow occurs in the scaled filter. 

It can be shown 19 that, in an addition operation, if two's comple- 
ment arithmetic is used, as is usually the case, then as long as the final 
result is within the representable numerical range, individual partial 
sums can be allowed to overflow without causing inaccuracies in the 
result. In this paper, it is assumed that all additions in a filter are 
done using two's complement arithmetic. Then, to guarantee that no 
distortion caused by overflow occurs at a cascade filter's output, only 
the input and output of each filter section need be constrained not to 
overflow. 

To simplify the discussion of scaling methods, the following defini- 
tions are used. Let 

Fi(z) = £ /<(*)*-*= UH } (z) (18) 

and 2 . . 1 ^ i ^ Ni 

Fi(z) = ifi(k)z- k = IT #,(*)• (19) 



;=i 



Also, let {Vi(n) } be the output sequence of F,(z) or ff<(»). Furthermore, 
assume that the maximum magnitude of numerical data representable 
in a filter is 1.0. Then the necessary overflow constraints on a cascade 
filter can be stated as 

|»<(n)|£ 1 1 ^ i£ N a , alln. (20) 

Necessary and sufficient conditions for (20) to hold for two classes 
of input signals are given below. Theorem 2 deals with the class of 
input sequences {x(n)\ which satisfies \x(n)\ ^ 1 for all n. For sim- 
plicity, this class is refererd to as "class 1." Theorem 3 deals with 
the class of inputs with transform X(e ju ) which satisfies 

1 f 2 * 

— / I X(e>") | da> ^ 1 . 
2tJ 

This class will be called "class 2." By virtue of the fact that 

1 f 2 " 
x ( n ) = — / X{e j «)e* n du (21) 

2ir J o 
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and hence 



\x(n)\ ^ — / \X(e'")\du, 

2t J o 



(22) 



class 2 is a subset of class 1. 



Theorem 2: Suppose \x(n)\ ^ 1. Then |y,-(n)| ^ 1, 1 ^ i ^ N, and 
all n, if and only if 



£ \Mk)\ = 1 i= 1, ■■■ ,N,. 



(23) 



Theorem 3: Suppose l/(2ir) f 2 ' | A'(e'») |dw ^ 1. Then \vi(n)\£ 1, 
1 ^ i ^ N, and all n, if and only if 



\Fi(e*»)\ = 1 i=l, 



, N, 



< u < 2tt. 



(24) 



Conditions (23) and (24) of Theorems 2 and 3 can be restated to 
give conditions on {£,-}. Recall that the i?,(z)'s are unique once H(z) 
is given, hence the F,-(z)'s and (/*(&) }'s are also unique. Equations 
(16), (18), and (19) give 



and 



/<(*) = [USj)Mk) 



Fiiz) =[US j )F i (z). 



(25) 
(26) 



Therefore, conditions (23) and (24) can be restated respectively as 

(27) 



and 



ri \s,\ ^ [l i/«(*)il 



II \S S \ ^ [ max lACe'-)!]" 1 



t = 1, ••• ,N,. 



(28) 



These then are conditions which, for the class of inputs concerned, a 
scaling method must satisfy. It will be shown next that in some sense 
optimum scaling methods are obtained when (27) and (28) are satis- 
fied with equality. For ease of reference, two scaling methods will 
first be defined. 

Define sum scaling to be the rule : 



r 2i 

n 8, = z 



/<(&) 



i = 1, ••• ,N, 



(29) 
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or stated recursively, 

.A-0 J 

Also, define peak scaling to be the rule: 



Si= ■ 



(30) 



i = 2, ■■• ,N.. 



II S f = [ max | A(e>") | ]-» t - 1, • • • , N. (31) 



or 



s«- , 



[ max | Pi(e ju ) | ] _I t - 1 



Kn«A max |A-(e'") 



i - 2, • • ■ , N, 



(32) 



Theorem 4: (riven an FIR transfer function to be realized in cascade 
form (as defined in Fig. 1) using fixed-point arithmetic of a given word- 
length, and given the ordering of filter sections, assume that: 

(i) The number of noise sources in each section (i.e., k { ) is inde- 
pendent of the scaling method, 
(ii) All filter coefficients can be represented to arbitrary precision. 
(Hi) No overflow is allowed to occur at the input and output of each 

section, 
(iv) The overall gain of the filter is maximized subject to no overflow 
at the filter output. 

Then each of the following scaling methods is optimum for the class of 
input signals staled, in the sense that it yields the minimum possible 
roundoff noise variance as defined in (13) among all scaling methods 
which satisfy conditions (Hi) and (iv) above for the class of inputs 
considered. 

(i) Sum scaling for class 1 signals, 
(ii) Peak scaling for class 2 signals. 

Thus optimal scaling methods are established for two classes of 
input signals. It is possible to define other classes of signals by con- 
sidering the "Lp norm" of their transforms. 1517 Specifically, the Lp 
norm of X(e ju ) is defined as 

IIWOL = [— r I W)| »d«l ' i ^ v ^ °° (33) 
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where for p = oo the limit as p — » °o of the right-hand side is meant. 
For each p, a class of signals can be defined consisting of those se- 
quences with transforms which satisfy 

ll*(«*OII,£ I- (34) 

Signals satisfying (34) will be referred to as Lp-norm constrained 
signals. Note that Li-norm constrained signals are simply class 2 
signals. 

For proofs of the following useful theorem, refer to Refs. 1, 20, and 21. 

Theorem 5: Let X(e iu ) and Y(e iu ) be transforms of sequences. Then 
(i) We''")!!. = max |X(«*»)| 

0Sa)g2ir 

(ii) l|*(*' w )W0||i ^ ||X(«*9UI *(«*•) ||. 

if l/p + \/q = 1, 1 ^ p,q ^«> 

(Hi) \\X(e^)\\ r ^ \\X(e»)\\, if l|r|«g». 

Since with input ( x (n) } , 

1 f 2 " 
Vi(n) = — / Fi(e' a )X(e^)e'^du, (35) 

2t Jo 

so that 

\vi(n)\ g — I \F<(e*)X(e*>)\dw = \\F t {e*>)X(e*)\\ h (36) 
Air J o 

by Theorem 5 (ii), 

\vi(n)\ ^\\Fi(e*)\\ p \\X(ei»)\\ q lgi^ N,. (37) 

Hence for L,-norm constrained signals, i.e., if ||X(e' u )|| g g 1, the 
following scaling method (L p -norm scaling) is obtained. 

Q 

11^.(^)11,-1 *"" Ar (38) 

i = 1, • • • , N„ 

or stated in terms of {Si}, 

II^=[||A(^)||p]-' i= 1, ■■■ ,N.. (39) 

Notice that by virtue of part (i) of Theorem 5, L M -norm scaling is 
just peak scaling which has been shown to be optimum for class 2, or 
Li-norm constrained, signals. Furthermore, by part (Hi) of the theorem, 
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the following hierarchy of classes of signals is obtained : 

class 1 D class 2 D L p -norm constrained D L g -norm 

constrained 

if 1 ^ V ^ Q ^ °° • 

In general, class 1 and class 2 signals are the most useful to consider. 
L 2 -norm constrained signals with L 2 -norm scaling are useful when all 
inputs to a filter have finite energy bounded by a known value. For, by 
Parseval's Theorem, the energy of {x(n)\ is simply given by 
(||X(e*»)|| 2 )*. Hence, if the input signals are first scaled so that their 
maximum energy is 1.0 (or the squared dynamic range of the filter), 
then Z,2-norm scaling is sufficient to ensure no overflow. 

Z/2-norm scaling finds greater application for FIR filters than for 
IIR filters because, in the former case, it is applicable for a larger class 
of input signals. In particular, an iVth-order FIR filter has only N 
samples of memory; thus if the input signal to an Nth-order FIR 
filter consists of bursts of energy spaced more than N samples apart 
with zero energy in between, then the filter will effectively "see" only 
one burst at a time. Hence, the maximum energy of a burst can be 
used as the bound on the energy of the input as far as scaling is con- 
cerned. Thus an infinite-energy signal can have the effect of a finite- 
energy signal on an FIR filter. 

Clearly, sum scaling and peak scaling can also be applied to IIR 
filters. 15 In fact, Theorems 2 and 3 are also valid for IIR filters. 
However, the input sequence needed in Theorem 2 to prove necessity 
in the case of IIR filters is an infinite-duration sequence extending to 
— oo with full dynamic range magnitudes, and signs that match those 
of [fi(k) } for some i. Since {/,(&) } for IIR filters is infinite in duration 
for all i, clearly such an input sequence is highly improbable. Hence, 
class 1 signals have been deemed too restrictive a description for 
ordinary inputs to an IIR filter, resulting in too stringent a scaling 

method. 15 

However, for FIR filters it is not difficult to find an input sequence 
within dynamic range which will require sum scaling to ensure no 
overflow, since only a small, finite portion of the sequence need match 
up with the j/,(/c)}'s. For example, if F^z) has a zero with angle w , 
tt/2 g wo < t, then all three samples of [fi(k)\ have the same sign; 
hence an input sequence need only have three consecutive samples of 
value 1 before |y x (n) | = L*|/i(fc) I for some n - 
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4.3 Properties of the Linear Phase Cascade Structure 

Two theorems regarding certain properties of the linear phase 
cascade form are now given. These results are useful in the investiga- 
tion of ordering of cascade filter sections. 14 

Theorem 6: Let Hi(z) be the transfer function for the ith section of a 
linear phase FIR filter in cascade form, where 

Hi(z) = b Qi + buz' 1 + bur*, 

and let co, be the angle of one of its zeros, —v ^ co, ^ tt. Then for all i: 



(?) max \Hi(e ju )\ = < 



\Hi(e*)\ Og | co, | <- 

2 



\Hi(e»)\ -^ | co, | g it 

2 

(ii) £ I bu | = max | #,(e>'<") | = max ( | #,(e'°) | , | H <(«'*) | ) . 



2 



1-0 



The next theorem is concerned with the equivalence of certain order- 
ings with regard to output noise variance. In particular, it states that 
with peak scaling each pair of sections in a filter which synthesize 
reciprocal zeros is completely interchangeable without affecting the 
output noise variance of the filter. With sum scaling, however, this 
is not necessarily true. Nevertheless, a weaker condition can be stated 
which says that, with sum scaling, if every pair of sections which 
synthesize reciprocal zeros of a filter is interchanged in position, then 
output noise variance is not changed. Figure 3 illustrates two such 
equivalent orderings. 

Theorem 7: Let |#,(z)| and {#, : (z)| be the section transfer functions of 
two orderings for a linear phase filter H (z) , both scaled by the same method, 
thus 

Hit) = ri#,(2) = n #;•(*). 

i=i .=1 

Then filters icith section transfer functions ji/,(z)} and (#,'(z)} produce 
identical output noise variances if either of the following conditions is 
true: 
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Fig. 3 — Two orderings with equal output noise variances. 

(i) Peak scaling is used and for each i, f/,(z) and flj(*) have either 
the same zeros or reciprocal zeros [i.e., Hi(zi) = OifHifa ) = Oj. 



(ii) Sum scaling is used and for all i, z { 
Zi is a zero of Hi(z). 



is a zero of Hi(z) whenever 
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